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This  paper  presents  theoretical  and  experimental  investigations  of  the  winter  operation  mode  of  a 
thermoelectric  cooling  and  heating  system  driven  by  a  heat  pipe  photovoltaic/thermal  (PV/T)  panel. 
And  the  energy  and  exergy  analysis  of  this  system  in  summer  and  winter  operation  modes  are  also  done. 
The  winter  operation  mode  of  this  system  is  tested  in  an  experimental  room  which  temperature  is 
controlled  at  18  °C.  The  results  indicate  the  average  coefficient  of  performance  (COP)  of  thermoelectric 
module  of  this  system  can  be  about  1.7,  the  electrical  efficiency  of  the  PV/T  panel  can  reach  16.7%,  and 
the  thermal  efficiency  of  this  system  can  reach  23.5%.  The  energy  and  exergy  analysis  show  the  energetic 
efficiency  of  the  system  in  summer  operation  mode  is  higher  than  that  of  it  in  winter  operation  mode,  but 
the  exergetic  efficiency  in  summer  operation  mode  is  lower  than  that  in  winter  operation  mode,  on  the 
contrary. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

In  the  past  decades,  because  of  the  development  of  the 
economic  and  industrial,  the  non-renewable  and  polluting  fossil 
fuels  were  consumed  largely.  And  it  is  necessary  to  explore  new 
energy  to  remit  the  needing  pressure  of  energy.  The  solar  energy 
is  a  new  way  to  relieve  the  pressure,  because  it  is  pollution  free 
and  easy  to  gain.  Many  devices  have  been  designed  to  gain  solar 
energy.  Flat-plate  solar  collector  is  steady  and  efficient  and  it  has 
been  widely  used  in  residential  water,  space  heating  and  commer¬ 
cial  or  industrial  applications.  Its  efficiency  can  be  improved  by 
reducing  its  size  and  obtaining  higher  temperature  fluid  at  outlet. 
There  are  many  different  highly-effective  techniques  which  have 
been  used  in  the  past  to  enhance  the  thermal  performance  of  solar 
collectors  including  the  methods  of  reducing  the  heat  loss  from  the 
top  surface  [1,2]  or  increasing  the  energy  gain  inside  the  solar  con¬ 
verter  [3,4],  The  building-integrated  dual-function  solar  collector  is 
a  new  structure  collector  which  can  perform  in  two  different 
modes:  working  as  a  passive  space  heating  collector  in  cold  sunny 
days  such  as  in  winter  or  working  as  a  facade  water  heating  collec¬ 
tor  in  hot  days  such  as  in  summer  [5-7], 
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The  PV/T  system  has  been  researched  widely  because  it  can  use 
solar  energy  more  sufficiently  to  produce  electrical  and  thermal 
energy  simultaneously.  Henning  Helmers  and  Korbinian  Kramer 
presented  a  performance  model  that  enables  yield  predictions  of 
hybrid  photovoltaic  and  thermal  (PVT)  collectors.  It  applied  for 
both  non-concentrating  (PVT)  and  concentrating  (CPVT)  systems. 
The  model  was  based  on  considerations  of  energy  balance,  heat 
transfer  and  the  dependence  of  the  photovoltaic  efficiency  on 
absorber  temperature  and  applied  to  measurement  data  of  a  CPVT 
collector  to  exemplify  the  procedure  and  to  validate  the  model  [8], 
Faizal  reported  the  energy,  economic  and  environmental  analysis 
of  metal  oxides  nanofluid  for  flat-plate  solar  collector  which  used 
nanofluid  as  working  fluid.  From  the  study,  it  was  estimated  that 
a  large  number  of  solar  collectors  can  be  saved  for  CuO,  Si02, 
Ti02  and  A1203  nanofluid.  The  average  value  of  220  MJ  embodied 
energy  can  be  saved  for  each  collector,  2.4  years  payback  period 
can  be  achieved  and  around  170  kg  less  C02  emissions  in  average 
can  be  offset  for  the  nanofluid  based  solar  collector  compared  to 
a  conventional  solar  collector  [9],  da  Silva  and  Fernandes 
researched  the  thermodynamic  modeling  of  hybrid  photovoltaic- 
thermal  (PV/T)  solar  systems,  pursuing  a  modular  strategy 
approach  provided  by  Simulink/Matlab.  And  the  results  showed 
that  the  modular  approach  strategy  provided  by  Matlab/Simulink 
was  applicable  to  solar  systems  modeling,  providing  better  code 
scalability,  faster  developing  time,  and  simpler  integration  with 


W.  He  et  al./Energy  Conversion  and  Management  84  (2014)  41-49 


Nomenclature 

A 

contact  area,  m2 

conductivity  of  the  thermoelectric  materials,  s/m 

Cp 

heat  capacit,  J/(kg  K) 

T 

Interval  time  of  data  collection,  s 

Epv 

Ete 

photovoltaic  output  power  per  unit  cell  area,  W/m2 
thermal  output  power  per  unit  collector  area,  W/m2 

{ 

PV  cell  packing,  factor 

Elrv 

photovoltaic  exergy  output  per  unit  PV  cell  area,  W/m2 

Subscripts 

E/tc 

thermal  exergy  output  per  unit  collector  area,  W/m2 

1 

ceramic  shell  of  thermoelectric  device 

EXsun 

exergy  input  of  solar  radiation,  W/m2 

2 

radiator 

gt 

solar  irradiance,  W/m2 

3 

foam  wall 

I 

current,  A 

4 

brick  wall 

K 

thermal  conductivity,  W/(kg  K) 

5 

copper  plate 

L 

thickness,  m 

6 

water  block 

M 

quality,  kg 

a 

ambient  environment 

P 

power,  W 

c 

cold  side  of  thermoelectric  module 

Q. 

Heat  flux,  W 

h 

hot  side  of  thermoelectric  module 

R 

thermal  resistance,  K/W 

heat-hr 

hot  side  of  thermoelectric  module  and  experimen¬ 

R' 

thermal  contact  resistance,  K/W 

tal  room 

R" 

contact  resistance  coefficient,  m2  K/W 

heat-cw 

cold  side  of  thermoelectric  module  and  heating  water 

S 

area,  m2 

n 

number 

T 

temperature,  K 

pipe-t 

top  of  heat  pipe 

AT 

temperature  variation  of  water,  °C 

pipe-m 

middle  of  heat  pipe 

U 

voltage,  V 

bottom  of  heat  pipe 

a 

absorptivity  or  Seebeck  coefficient  thermoelectric 

panel 

absorber  panel 

materials,  V/K 

experimental  room 

fl 

Boltzmann’s  constant,  W/(m2  K4) 

s 

thermoelectric  module 

£ 

Exergetic  efficiency 

solar 

PV  cell 

the 

thermal  efficiency  of  the  system 

[ 

thermoelement 

t\pv 

electrical  efficiency  of  solar  cells 

te 

thermal  energy 

tjpvt 

total  efficiency  of  the  system 
emissivity  of  absorber  panel 

w 

water  in  the  storage  tank 

external  computational  tools,  when  compared  with  traditional 
imperative-oriented  programming  languages  [10],  Highly  thermal 
efficiency  can  be  achieved  via  the  combination  of  solar  collector 
and  other  devices.  The  composite  structure  of  solar  collector  and 
heat  pump  is  common  one.  Jose  Fernandez-Seara  and  Carolina 
Pineiro  reported  the  experimental  analysis  of  a  direct  expansion 
solar  assisted  heat  pump  with  integral  storage  tank  for  domestic 
water  heating  under  zero  solar  radiation  conditions  [11].  Moren- 
o-Rodriguez  and  Gonzalez-Gil  presented  the  theoretical  model 
and  experimental  validation  of  a  direct-expansion  solar  assisted 
heat  pump  for  domestic  hot  water  (DHW)  applications.  The 
acquired  experimental  coefficient  of  performance  was  found  to 
be  in  the  rank  of  1.7-2.9.The  DHW  tank  temperature  over  the 
course  of  the  study  is  51  °C[12],  S.K.  Chaturvedi  researched  the 
solar-assisted  heat  pump  which  was  sustainable  for  low-tempera¬ 
ture  water  heating  applications.  Results  indicated  that  the  DX- 
SAHP  (Direct  expansion  solar  assisted  heat  pump)  water  heaters 
systems  when  compared  to  the  conventional  electrical  water  heat¬ 
ers  were  both  economical  as  well  as  energy  conserving.  The  analy¬ 
sis  also  revealed  that  the  minimum  value  of  the  system  life  cycle 
cost  was  achieved  at  optimal  values  of  the  solar  collector  area  as 
well  as  the  compressor  displacement  capacity  [13],  Gang  and 
Huide  presented  a  dynamic  model  of  a  heat  pipe  PV/T  system 
and  constructed  a  test  rig.  Experiments  were  conducted  to  validate 
the  results  of  the  simulation.  Based  on  the  validated  module,  the 
performances  of  the  heat  pipe  PV/T  system  were  studied  under  dif¬ 
ferent  parametric  conditions,  such  as  water  flow  rates,  PV  cell  cov¬ 
ering  factor  of  the  collector,  tube  space  of  heat  pipes,  and  kinds  of 
solar  absorptive  coatings  of  the  absorber  plate  [14], 

Solar  thermoelectric  refrigerator  were  reported  by  many 
researchers  [15-17],  but  thermoelectric  heating  was  not  a  common 
research  field.  This  article  introduced  the  application  of  combining 


of  PV/T  system,  thermoelectric  modules  and  building.  The  summer 
operation  mode  of  the  system  which  can  provide  cooling  for  room 
and  heat  domestic  hot  water  for  user  was  tested  with  the  small- 
scale  version  [18],  The  results  indicate  the  system  has  a  higher 
electrical  efficiency  (15.4%)  and  thermal  efficiency  (29.4%),  and 
the  thermoelectric  module  has  a  strong  cooling  ability  that  COP 
is  about  0.45.  Now  the  experimental  performance  and  theoretical 
analysis  of  the  system  in  winter  operation  mode  are  presented  in 
this  paper.  And  energy  and  exergy  analysis  are  also  done  in  sum¬ 
mer  and  winter  operation  modes,  which  are  validated  by  the 
small-scale  system  and  this  paper. 

2.  Experimental  structure  and  working  principles  of  the  system 

2.1.  Experimental  structure 

The  experimental  structure  of  the  thermoelectric  heating  sys¬ 
tem  driven  by  a  heat  pipe  PV/T  panel  is  shown  in  Figs.  1  and  2. 

The  system  is  composed  of  a  heat  pipe  PV/T  panel,  model  of 
maximum  power  point  tracking  and  controling  the  solar  charge, 
experimental  room,  thermoelectric  modules,  heat  exchangers, 
pumps,  fans  and  storage  tank.  The  experiment  was  carried  out  on 
the  sunny  days  of  winter.  The  dates  which  were  recorded  include 
the  temperature  of  the  PV/T  panel,  heat  pipe,  the  cold  side  and 
hot  side  of  thermoelectric  modules,  the  experimental  room  and 
comparative  room,  ambient  and  the  water  of  storage  tank,  the  solar 
irradiance,  the  output  Voltage  and  current  of  the  PV/T  system,  the 
input  Voltages  and  currents  of  thermoelectric  modules,  the 
amount  electricity  used  in  the  experimental  and  comparative  room 
in  which  EHS  used.  The  difference  of  the  mathematical  models 
between  the  summer  operation  mode  and  winter  operation  mode 
is  that  the  direction  of  current  input  is  opposite,  hence  the  hot  side 
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Fig.  2.  Experimental  installations. 


in  summer  which  release  heat  to  cooling  water  is  the  cold  side  in 
winter  which  absorbs  heat  from  the  heating  water,  and  the  cold 
side  in  summer  which  absorbs  heat  from  the  experimental  room 
is  the  hot  side  in  winter  which  release  heat  to  room.  The  difference 
of  the  system  in  winter  operation  mode  is  that  there  are  four  ther¬ 
moelectric  modules  (six  modules  in  summer)  and  30  L  (150  L  in 
summer)  water  is  used  in  the  experiment,  and  the  temperature 
of  experimental  in  winter  is  set  at  18  °C  (26  °C  in  summer). 


2.1.1.  Heat  pipe  PV/T  panel 

The  effective  irradiation-collection  area  of  the  HP-PV/T  solar 
collector  is  1.95  m2,  and  the  PV  cell  area  is  1.53  m2.  The  structure 


Table  1 

The  performance  parameters  of  it. 

Heat  pipe  PV/T  panel  ^  Pmwl  (W)  Ump  (V)  Imp  (A)  U0C(V)  fsc(A) 

18%  272  49.5  5.49  59.4  5.94 


of  it  is  showed  in  Fig.  3.  The  PV  modules  comprise  98  pieces  of 
small  PV  cells  that  are  arranged  in  series  with  dimensions 
ofl25  x  125  mm,  and  its  photovoltaic  characteristics  under  stan¬ 
dard  test  conditions  (at  irradiationintensity  of  1000  W/m2  and 
temperature  of  25  °C)are  presented  in  Table  1. 

2.1.2.  Model  of  maximum  power  point  tracking  solar  charge  controller 

The  model  of  maximum  power  point  tracking  solar  charge  con¬ 
troller  (MPPT)  is  Tracer-221 0RN.  The  performance  parameters  of  it 
are  showed  in  Table  2. 

The  tracer  series  controller  is  for  off-grid  solar  systems  and  it 
controls  the  charging  and  discharging  of  the  batteiy.  It  is  especially 
suitable  for  the  thermoelectric  modules  systems.  The  controller 
features  a  smart  tracking  algorithm  inside  that  maximizes  the 
energy  from  the  solar  PV  module  and  charges  the  battery.  The 
low  voltage  disconnect  function  prevents  the  battery  from  over 
discharging. 

2.1.3.  Indoor  temperature  control  system  within  the  experimental  and 
comparative  rooms 

The  indoor  temperature  controller  system  in  the  experimental 
room  and  comparative  room  consist  of  two  parts:  an 
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Performance  parameters  of  MPPT. 


MPPT  System  voltage 

1 2/24VDC20A20A 


Rated  battery  current  rate 
100VDC 


Max.  PV  input  voltage  Max.  PV  input  Power 


air-conditioner  system  that  supplies  cooling  capacity  and  records 
the  amount  of  heat  removed  from  the  indoor  area  and  an  electrical 
heating  system  (EHS)  that  supplies  heating  capacity  and  records 
the  amount  of  heat  released  into  the  indoor  area.  These  two  sys¬ 
tems  work  together  to  maintain  the  temperature  at  desired  setting. 
The  heating  ability  of  the  thermoelectric  modules  can  be  calculated 
by  the  difference  between  the  heats  removed  and  released. 


lectric  module 


2.2.  Working  principle  of  the  system  in  winter 


Fig.  4.  Construct  of  the  cooling  system. 


The  working  principle  of  the  system  in  winter  is  showed  in  Fig.  1 . 
The  heat  pipe  PV/T  panel  (1 )  is  located  outside  of  the  experimental 
room.  It  is  combined  with  the  heat  exchanger  (2)  via  heat  pipes. 
Inside  of  the  experimental  room,  the  hot  sides  of  the  thermoelectric 
modules  (5)  are  connected  to  the  radiators  (6),  and  the  cold  sides  of 
the  modules  (7)are  connected  with  the  water  blocks  (8).The  water 
blocks  (8),  storage  tank  (4)  and  heat  exchanger  (2)are  connected 
together  via  circuit  (3).  When  the  system  works,  the  water  flows 
through  heat  exchanger  (2),  water  blocks  (8)  and  then  back  to  stor¬ 
age  tank  (4).  The  heat  pipe  PV/T  panel  (1)  heats  the  circulating 
water  in  the  heat  exchanger  (2)  via  heat  pipes,  and  then  the  cold 
sides  of  thermoelectric  modules  (5)  absorb  heat  from  the  circulat¬ 
ing  water.  The  hot  sides  of  thermoelectric  modules  (5)  release  heat 
to  the  experimental  room  via  radiators  (6)  and  fans. 

3.  Mathematical  models  and  error  analysis 


Here,  Qh  is  the  heat  released  by  the  hot  side  of  thermoelectric  mod¬ 
ule.  lift  is  the  input  power  of  the  thermoelectric  module. 

3.4.  Electrical  efficiency  of  the  heat  pipe  PV/T  [19] 

%v  =  »7rro/  [1  -  0.0045(Tso/ar  -  298.15)]  (4) 

Here,  r]pv  is  the  electrical  efficiency.  rjT ^  is  the  module’s  electrical 
efficiency  at  the  reference  temperature  which  is  0.18.  Tsolar  is  the 
module  operating  temperature. 

3.5.  The  thermal  efficiency  of  the  system 


Qtan*  +  S£^&»i 

T^GnTSpan 


(5) 


3.1.  The  heat  balance  between  the  hot  side  of  thermoelectric  module 
and  experimental  room  is  expressed  as  follows 


Q-heat-hr 


Th~Tr 
Rj  +  R'u  +  R2 


0) 


Here,  Qheat-hr  is  the  heat  transfer  between  the  hot  side  of  the 
thermoelectric  module  and  the  experimental  room,  and  Th  is 
the  temperature  of  the  hot  side  of  thermoelectric  module,  Tr  is 
the  temperature  of  the  experimental  room.  R,  is  the  thermal  resis¬ 
tance  of  the  ceramic  shell  of  thermoelectric  module.  R'12  is  the  ther¬ 
mal  contact  resistance  between  the  thermoelectric  module  and 
radiator,  R2  is  the  thermal  resistance  of  the  radiator. 


3.2.  The  construct  of  the  cooling  system  is  showed  in  Fig.  4 


Here,  t]te  is  the  thermal  efficiency  of  the  system.  Qu,nk  is  the  heat 
that  the  water  in  storage  tank  has  gained.  XS”tornsQ.w  is  the  heat 
which  the  thermoelectric  modules  released.  Qtank  +  Jf,nZ"s°u,rnsQ.hi  is 
the  total  heat  which  the  system  has  released.  CTi  is  the  solar  radia¬ 
tion  flux  (irradiance)  on  module  plane,  t  is  the  time  interval  of  data 
collection.  Spanel  is  the  area  of  the  absorber  plate.  nsolar  is  the  number 
of  the  data  has  been  recorded.  X)";?”'“rGTiTSp0n!,/  is  the  total  energy 
which  the  PV/T  panel  has  gained. 

3.6.  The  energy  and  exergy  analysis  of  the  heat  pipe  PV/T  system 

From  the  first  law  of  thermodynamics,  the  energetic  efficiency 
of  the  system  r\pvt  is  composed  by  two  parts:  the  electrical  effi¬ 
ciency  ripv  and  thermal  efficiency  ijte,  and  the  equation  are  given 
as  follows  [20]: 


The  heat  balance  between  the  cold  side  of  thermoelectric  mod¬ 
ule  and  heating  water  is  expressed  as  follows: 

Q'"Mt~cw  =  Ri  +  R5  +  Re  +  R\  s  +  R'56  ^ 

Here,  Qheat-cw  is  the  heat  transfer  between  the  cold  side  of  thermo¬ 
electric  module  and  heating  water.  Tw  is  the  temperature  of  heating 
water.  Tc  is  the  temperature  of  the  cold  side  of  thermoelectric  mod¬ 
ule,  R5  is  the  thermal  resistance  of  copper  plate.  R6  is  the  thermal 
resistance  of  wall  of  water  block.  R'15  is  the  thermal  contact  resis¬ 
tance  between  the  ceramic  and  copper  plate.  R'56  is  the  thermal  con¬ 
tact  resistance  between  the  copper  plate  and  water  block  (see  Fig.  4). 

3.3.  Coefficient  of  performance  (COP)  of  the  thermoelectric  module 
C0P  =  W-t  (3) 


ft*  (SpanelEte  +  SsoiarEpv)dt 
==  SpanelftGrdt 


Z%v  +  Vte 


(6) 


(7) 


Here,  Epv  is  the  photovoltaic  output  power  per  unit  cell  area.  Ea  is 
the  thermal  output  power  per  unit  collector  area.  Gris  the  solar  irra¬ 
diation  per  unit  collector  area.  {  is  the  PV  cell  packing  factor. 

The  first  law  of  thermodynamics  assessed  the  performance  of 
the  PV/T  system,  however,  these  assessments  cannot  perfectly 
describe  the  performance  of  it,  because  only  under  the  condition 
of  existing  a  temperature  difference  between  a  high-temperature 
heat  source  and  a  low-temperature  heat  sink,  can  the  thermal 
energy  produce  work.  Thus,  the  exergy  analysis  was  based  on  the 
second  law  of  thermodynamics,  which  revealed  a  system  with  a 
reasonable  degree  of  energy  and  can  better  evaluate  the  perfor- 
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mance  of  the  PV/T  system.  The  exergetic  efficiency  can  be  defined 


as  follows  [21]: 

fi(SPanelEX*+SsolarE(ppv)dt  ^  „ 

SpanelftEXsundt  ~~ m  £te 

(8) 

EXpv  =  Epv 

(9) 

(10) 

(11) 

Here,  Eypv  is  the  photovoltaic  exergy  output  per  unit  PV  cell  area. 
Eyle  is  the  thermal  exergy  output  per  unit  collector  area.  EXsun  is 
the  exergy  input  of  solar  radiation.  Tsun  is  the  solar  radiation  tem¬ 
perature  which  is  6000  K. 


3.7.  The  energy  transfer  process  of  the  system 


The  current  of  the  thermoelectric  module  [22]: 

„  f  In  In  \  2 ntlt 

t~nt[(E„Sn  +  apSp)  ~tr„Sn 


(12) 


_U,  _  U,(TnSn 
‘  ~  Rt  ~  2  ntLt 


(13) 


Here,  Qt  is  the  resistance  of  the  thermoelectric  module,  a  is  the 
Seebeck  coefficient  of  thermoelectric  materials  (the  p-type  and 
the  n-type  assumed  to  be  same),  It  is  the  current,  St  is  the  cross- 
sectional  area,  It  is  length  of  thermoelement,  a  is  the  conductivity 
of  the  thermoelectric  materials,  nt  is  the  number  of  thermoelements 
Ut  is  the  voltage  of  thermoelectric  module. 

The  temperature  of  hot  side  of  thermoelectric  module: 

Qw  A  7V  '=*  Th 2  =  Tr  +  QhlRhot  (14) 


Here,  Rhot  is  the  total  thermal  resistance  between  the  hot  side  of 
thermoelectric  module  and  experimental  room. 

The  temperature  of  cold  side  of  thermoelectric  module: 

del  =Jw'R  Ta  =>  Ta  =  TW1  -  QclRcold  (15) 


Here,  RcoW  is  the  total  thermal  resistance  between  the  cold  side  of 
thermoelectric  module  and  water  block. 

The  temperature  of  water  in  storage  tank: 


■t(upumpIpumpr\pump  +  npipeT-^E  _  nsQc  -  Q,oss) 


Here,  t  is  the  interval  time  of  data  collection.  Upumpipumpt]pump  is  the 
heat  which  the  pump  released.  Rpw  is  the  total  thermal  resistance 
between  the  heat  pipe  PV/T  panel  and  water  in  the  heat  exchanger. 
npipe  is  the  heat  which  was  transferred  from  PV/T  panel  to 

the  water  in  storage  tank.  nsQc  is  the  heat  which  the  thermoelec¬ 
tric  modules  absorbed.  Q,ms  is  the  total  heat  loss  of  storage  tank 
and  circulation  line.  t(Upum!,Ipumptjpump  +  nm  Tp“"PpJ'“'  -  nsQc  -  Qlo„) 
is  the  heat  which  storage  tank  gained  in  the  interval  time  of  data 
collection.  A  flow  chart  of  this  system  is  showed  in  Fig.  5. 

According  to  the  first  section  in  Fig.  6,  with  the  voltage  (Upvt) 
and  current  (fpvt)  which  the  heat  pipe  PV/T  panel  produces,  actual 
power  (Pen p)  of  PV/T  system  can  be  calculated,  and  then  it  will  be 
compared  with  the  theoretical  power  ( Psijn )  which  can  be  calcu¬ 


lated  with  Eq.  (4).  If  the  difference  is  within  the  setting  range,  it 
will  continue  to  the  second  section,  else  it  will  return  to  the 
beginning. 

During  the  second  section,  the  current  (In)  of  the  thermoelectric 
modules  can  be  calculated  with  the  voltage  (Utl)  (Eqs.  (12)  and 
(13))  which  the  battery  output.  With  the  Tciexp  and  Thlexp,  the  heat 
absorbed  by  cold  side  (Lislm  and  heat  released  by  hot  side  Qhlsim  of 
thermoelectric  module  can  be  calculated.  With  the  temperature 
of  experimental  room  Tr  and  dhi^,  the  Th2pm  can  be  calculated 
(equation  14).  With  the  temperature  of  water  Twlexp,  and  the 
QcW  Q-heat  pwi exp ■  Tc2tlm  and  Tw2sim  be  calculated  (Eqs.  (15)  and 
(16)),  then  Tc2sim,  Th2slm  and  Tw2jim  will  be  compared  with  the  exper¬ 
imental  values.  If  the  difference  is  within  the  setting  range,  the 
program  will  be  finished,  else  it  will  return  to  the  beginning. 


3.8.  Error  analysis 

The  experimental  error  of  the  independent  variables,  such  as 
temperature,  output  electricity  and  solar  irradiation,  is  determined 
by  the  accuracy  of  the  corresponding  instrument.  While  the  exper¬ 
imental  error  of  the  dependent  variables,  including  the  overall  sys¬ 
tem  heat  gain,  r/pvt,  r\pv  and  t/te  can  be  calculated  from  the 
experimental  error  of  the  independent  variables  according  to  the 
theory  of  error  propagation. 

The  relative  error  (RE)  of  the  dependent  variable  y  is  calculated 
as  follows  [21]: 

RE_d +  (1„ 

y  ax,  y  dx2  y  x„  y 


y  =  J  (xi,x2,...x„) 


(18) 


where  xf,  (i  =  1  ...n)  is  the  variable  of  the  dependent  variable  y,  and 
0//8x  is  the  error  transferring  coefficient  of  the  variables. 

The  experimental  relative  mean  error  (RME)  during  the  test  per¬ 
iod  can  be  expressed  as: 


E?l*cl 


(19) 


Based  on  the  Eqs.  (21)-(23),  The  RME  of  all  variables  discussed  is 
listed  in  Table  3. 


4.  Results  and  analysis 

In  the  morning,  with  the  solar  irradiance  increasing  from 
200  w/m2  to  700  w/m2,  the  output  power  of  the  heat  pipe  PV/T 
panel  increased  fast  from  50  w  to  170  w  (shown  in  Fig.  6),  and 
the  temperature  of  the  PV/T  panel  increases  steady,  hence  the  heat 
which  the  water  has  absorbed  and  the  temperature  of  it  both 
increase.  The  temperature  of  cold  side  of  the  thermoelectric  mod¬ 
ule  as  a  similar  trend  with  the  water,  because  it  absorbs  heat  from 
the  water,  but  the  temperature  of  hot  side  of  the  thermoelectric 
module  does  not  change  obviously,  since  the  heat  which  the  cold 
side  absorbs  (As  shown  in  Fig.  8,  the  difference  of  temperature 
between  the  water  and  the  cold  side  is  very  steady,  hence  the  heat 
transferred  between  them  do  not  change  nearly)  and  the  input 
power  of  the  thermoelectric  module  which  has  a  steady  input  volt¬ 
age  supplied  by  battery  (shown  in  Fig.  9)  are  both  steady,  and 
relationship  between  the  heat  release  and  the  temperature  of  the 
hot  side  is  nearly  linear  correlation  (the  current  is  steady).ln  the 
afternoon,  with  the  solar  irradiance  decreasing  from  700  w/m2  to 
300  w/m2,  the  output  power  and  the  temperature  of  PV/T  panel 
both  decrease,  but  the  temperature  of  panel  (average  35  °C)is 
always  higher  than  the  water  (average  21  °C),  hence  the  tempera¬ 
ture  of  water  and  cold  side  of  thermoelectric  module  keep  increas- 
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Fig.  5.  Flow  chart  of  simulation. 


ing  slowly,  and  the  temperature  of  hot  side  keeps  steady.  The  COP 
is  stable  before  11:00,  and  the  reason  is  that  the  increasing  rate  of 
Qh  is  similar  with  input  power  of  thermoelectric  module.  With  the 


temperature  of  water  increasing,  Tc  and  Qc  also  increase,  and  the 
rate  of  increase  is  much  faster.  Qj,  contains  two  parts:  the  Qc  and 
UtIt.  The  UJt  is  very  steady  from  11:00  to  13:00,  and  then  begins 
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Fig.  6.  Variations  of  output  power  and  solar  irradiance  on  a  sunny  day. 


Fig.  9.  Variations  in  the  input  power,  COP  and  Q,,  of  the  thermoelectric  module. 


Table  3 

Experimental  RME  of  the  variables. 

Variable  T  CT  ;/„  i/pv  ;/pvt 

RME  0.0693%  1%  0.4%  0.09%  1.49% 


*  T(panel) 

*  T(pipe-M) 
°  T(pipe-T) 

*  T(pipe-B) 


Fig.  7.  Variations  of  the  temperature  of  PV/T  panel  and  heat  pipe  on  a  sunny  day. 


8:00  9:00  10:00  11:00  12:00  13:00  14:00  15:00  16:00 

Time 

Fig.  8.  Variations  in  Th,  Tc,  T„  and  ambient  temperature  Tamb  on  a  sunny  day. 


to  decrease.  Qc  increases  steadily  from  11:00  to  16:00.  Hence  the 
COP  (COP  =  §f  —  ‘'y  )  begins  to  increase  from  11:00. 

Fig.  7  shows  the  variation  in  temperature  of  heat  pipe  PV/T 
panel  and  heat  pipe  through  all  day.  The  temperature  of  the  heat 
pipes  in  the  bottom  is  highest,  because  it  is  the  evaporator  which 
absorbs  heat  from  the  panel.  The  temperature  of  the  heat  pipe  at 
the  top  is  lowest,  because  it  is  the  condenser  which  releases  heat 
into  the  water.  The  temperature  in  the  middle  is  between  that  in 
the  bottom  and  at  the  top,  because  the  vapor  turns  into  water 
when  its  heat  is  absorbed  by  the  circulating  water,  and  then  the 
water  returns  back  to  the  bottom  of  heat  pipe  with  the  force  of 
gravity  which  will  decrease  the  temperature  of  middle  during  the 
process. 

As  shown  in  The  Fig.  10,  the  heat  released  by  electrical  heating 
system  in  experimental  room  and  comparative  room  can  be  calcu¬ 
lated.  And  all  the  data  is  shown  in  Table  4. 

The  simulated  and  experimental  electrical  efficiency  of  the  heat 
pipe  PV/T  panel  are  showed  in  Fig.  11.  The  electrical  efficiency  is 
related  to  the  temperature  of  the  panel.  The  temperature  of  panel 
increases  in  the  morning,  peaking  at  noon.  It  decreases  during  the 
afternoon  (Fig.  6).  Electrical  efficiency  decreases  in  the  morning, 
reaches  its  minimum  at  noon,  and  then  begins  to  rise  in  the  after¬ 
noon.  The  simulated  and  experimental  electrical  efficiency  are 
founded  to  be  correlated  with  each  other,  and  the  average  values 
arel7.6%  and  16.7%. 


Fig.  10.  Variation  in  the  amount  electricity  used  in  the  experimental  and 
comparative  room  in  which  EHS  used. 
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Heat  release  compariso 


Heat  released 
of  EHS 


Heat  transferred 
from  ambient 


Experiment  1.20(KWh)  0.735(KWh)  -1.935(KWh) 

Comparative  1.855  (KWh)  -1.855(KWh) 


Fig.  11.  Electrical  efficiency  of  the  heat  pipe  PV/T  panel. 


As  shown  in  Figs.  12  and  13,  all  the  simulative  and  experimental 
temperature  of  cold  side,  hot  side  of  thermoelectric  module  and 
water  in  storage  tank  are  corrected  with  each  other,  and  the  max¬ 
imum  deviation  is  lower  than  1  °C.  The  temperature  of  water  in 
storage  tank  increases  byl5.2°C,  and  the  simulated  and  experi¬ 
mental  thermal  efficiencyare23.9%and  23.5%,  respectively. 

The  thermal  efficiency  of  the  system  in  winter  operation  mode 
is  lower  than  it  in  summer  operation  mode.  The  main  reason 
includes  two  factors:  (1)  there  are  six  thermoelectric  modules 
(total  input  power  156W)  used  in  summer  which  are  powered 
by  the  heat  pipe  PV/T  system  (average  output  power  161  W).  There 
are  four  thermoelectric  modules  (total  input  power  104  W)  used  in 
winter,  and  the  average  output  power  of  the  heat  pipe  PV/T  system 
is  125  W.  There  is  some  electricity  which  is  not  used  and  stored  in 


Fig.  13.  Temperature  of  the  hot  side  of  thermoelectric  module. 

the  battery,  hence  the  heat  thermoelectric  modules  releases 
decreased:  (2)  the  thermal  contact  resistance  is  related  to  the  tem¬ 
perature,  and  while  the  temperature  is  increasing,  the  thermal  con¬ 
tact  resistance  is  decreasing  [23],  The  temperature  of  the 
circulating  water  in  winter  operation  mode( average  14  °C)  is  much 
lower  than  that  in  summer  operation  mode(average  31  °C),  hence 
the  ability  of  heat  transfer  between  the  heat  pipe  PV/T  panel  and 
the  circulating  water  in  winter  operation  mode  is  worse  than  that 
in  summer  operation  mode. 

The  energetic  efficiency  of  the  heat  pipe  PV/T  system  in  summer 
operation  mode  (41.49%)  is  higher  than  that  in  winter  operation 
mode  (36.6%),  but  the  exergetic  efficiency  of  the  heat  pipe  PV/T 
system  in  summer  operation  mode  (12.9%)  is  lower  than  that  in 
winter  operation  mode  (14.14%).The  main  influencing  factor  is 
the  ambition  temperature.  The  ambition  temperature  in  winter 
(9  °C)  is  lower  20°Cthan  that  in  summer  (29  °C),  hence  the  temper¬ 
ature  of  the  panel  in  winter  is  lower  than  that  in  summer,  and  elec¬ 
trical  efficiency  of  heat  pipe  PV/T  system  in  winter  (16.7%)  is  higher 
1.3%  than  that  in  summer  (15.4%).  The  thermal  exergetic  efficiency 
in  winter  (1.08%)  is  higher  0.27  than  that  in  summer  (0.81%). 

5.  The  energy  analysis  of  the  heat  pipe  PV/T  system  in  summer 
and  winter 

The  energy  analysis  presented  in  this  section  is  mainly  based  on 
the  first  law  of  thermodynamics.  And  the  theoretical  model  based 
on  thermal  energy  balance  is  employed  for  the  study  of  the  heat 
pipe  PV/T  system: 

Q-absorbed  =  Qaccumulated  +  Pused  +  Qjosr  (20) 

The  absorbed  energy: 

Q-absorbed  =  GrSpanel  (21) 

The  accumulated  energy: 

Qaccumulated  =  ^p,cMpanel  (22) 

Here,  Cp,c  and  Mpalw,  are  the  specific  heat  capacity  and  quality  of 

absorber  panel,  respectively. 

The  useful  energy: 

Pused  =  t7rre/SsofarGr[l  -  0.0045(Tpane,  -  298.15)]  +  nptpJpane^Tw 

(23) 


Fig.  12.  Temperature  of  the  cold  side  of  thermoelectric  module  and  the  water  in  the  Here,  S,„lar  is  the  area  of  solar  array  to  receive  solar  irradiation, 

storage  tank.  t/7raf  SSOtarGT[l  -  0.0045(Tpane/  -  298.15)]  is  the  energy  which  the  solar 
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and  heating  capacity  about  127.97  KWh  in  winter  season 
in  this  case. 

(3)  As  the  ambition  temperature  and  operation  temperature  are 
important  factors,  the  energetic  efficiency  of  this  system  in 
summer  operation  mode  is  higher  than  that  in  winter  oper¬ 
ation  mode,  but  the  exergetic  efficiency  of  the  heat  pipe  PV/T 
system  in  summer  operation  mode  is  lower  than  that  in  win¬ 
ter  operation  mode. 
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ergy.  npipe  Tpa^pJw'  is  the  heat  which 


Apane,  ~  Kane, )  +  (2-8  +  3.0 U)Spanel(Tpanel  - 


’amt) 

(24) 

(I  is  the  Boltzmann’s  constant.  Kpa„ei  is  the  emissivity  of  absorber 
panel.  KpaneiSpaneill(T^ariel  -  T^mb)  is  the  radiation  heat  transferred 
between  the  panel  and  ambition.  (2.8  +  3.0u)Spanei(Tpa„et  -  Tamt)  is 
the  convective  heat  transfer  between  the  panel  and  ambition. 

Buildings  composited  with  thermoelectric  cooling  and  heating 
systems  driven  by  a  heat  pipe  PV/T  system  can  reduce  cooling  load 
in  summer  operation  mode  and  heating  load  in  winter  operation 
mode.  The  electrical  efficiency,  total  cooling  capacity  in  summer 
and  heating  capacity  in  winter  are  calculated  with  the  climate 
design  data  2005  Hefei  of  china  (Climate  Design  Data  2005  ASHRAE 
Handbook),  and  they  are  showed  in  Fig.  14  and  Table  5. 


6.  Conclusion 


Theoretical  and  experimental  investigation  on  the  winter  oper¬ 
ation  mode  of  a  thermoelectric  cooling  and  heating  system  driven 
by  heat  pipe  photovoltaic/thermal  panel  are  presented  in  this 
paper.  And  the  energy  and  exergy  analysis  of  this  system  in  summer 
and  winter  operation  mode  are  also  done.  The  result  indicates  that: 

(1 )  In  this  case,  the  system  has  a  good  heating  capacity,  the  elec¬ 
trical  efficiency  of  the  heat  pipe  PV/T  panel  is  about  16.7%, 
the  thermal  efficiency  of  the  system  is  about  23.5%. 

(2)  Based  on  calculation  with  the  climate  design  data  2005  Hefei 
of  China,  one  building  composited  with  this  system  can  pro¬ 
vide  cooling  capacity  about  35.77  KW  h  in  summer  season 


